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I. INTRODUCTION 


The search for new quark flavors remains an important and interesting task. The ex¬ 
istence of a fourth generation cannot be excluded on the basis of present measurements. 
The p-parameter only restricts the possible mass-splitting and the limit on the ‘number of 
neutrinos’ refers to light ones only. In fact, some GUTs even require quarks that do not fit 
into the standard model scheme. 

Because charm as well as bottom have been discovered through hadronic production of 
their quarkonium bound states J/'tp and T respectively, some authors [|,|| examined the 
prospects for discovering new flavors at future hadron colliders through a similar mechanism. 
Given favorable circumstances, gluon fusion would be the dominant source for quarkonium 
production and therefore especially the pseudoscalar ground state rj could be produced with 
sufficient rate. Also the ^Si state -0 might be accessible. 

For the lighter member of a fourth generation doublet the single quark decay, i.e. the 
decay of one of the constituents in the quarkonium, is likely to be suppressed due to small 
intergeneration mixing angles. Thus one may in a hrst step ignore this channel and only 
consider the annihilation decays. The latter include some new and distinctive modes which 
may even become dominant. The most important one could he rj ^ ZH, which might even 
offer a way to discover the Higgs boson |Q,^. 

Predictions for the decays p,0 — 77 , 7 Z, ZZ, WW, ZH have been derived in Born ap¬ 
proximation in references . QGD corrections, however, are only partly known. The aim 
of this work is to £11 this gap. 

The paper is organized as follows: The calculational method employed can be found in 

and will be discussed only briefly in the following section II, together with some general 
considerations. In section III the results of our calculations will be presented. Gompact 
approximations will be given in section IV, a brief summary in section V. 


II. GENERAL CONSIDERATIONS 

Not yet discovered quarks must be heavy and a nonrelativistic treatment of their bound 
states is adequate. As is well known, in this case the decay width of an S-wave bound state 
factorizes into a nonperturbative part - the wavefunction at the origin - and a perturbative 
part which is proportional to the free quark scattering cross section: 

( 1 ) 

where v denotes the relative velocity of the qq system and M. the free scattering amplitude. 
The factorization in this form applies for S-states only (for P-states see [^) and to the order 
considered in this paper. Relativistic corrections first enter at (P(a|), not considered in this 
work. 

There are several ways to calculate the rate: one is to simply compute the spin averaged 
cross section (with a modified statistical factor 1/4 1/(2*S' -|- 1) where S refers to the spin 

of the bound state), if only the desired spin configurations can contribute to the sum. This 
may require care concerning possible couplings. Gonsider for example the yiif-mode: among 
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the S waves only the spin 1 state with = 1 can decay that way, and the approach is 
straightforward. However, if we switch to the ZhT-mode, both states with = O”"*" and 
= 1- are present in the sum, albeit with unknown relative weight. In this case one 
may identify the parity violating part oc a • f in the amplitude with the r] decay, the part 
oc + a^) with -0. A second way is to project the appropriate amplitudes using the method 
derived in |]I|. Both methods were used to obtain and check the results given below, with 
the exception of the decays into two hhs. There only the second one was applied because 
the separation of the various couplings is inconvenient. 

The calculation of the wavefunction Rs requires the knowledge of the QCD potential. 
To get rid of the dependence on the potential model all widths are normalized to 

rei,^77) = 12|fis(0)P^ (2) 

and only the ratios 


r(X ^ ah) 


( 3 ) 


are presented. 

The zeroth order generic Feynman diagrams responsible for the annihilation decays are 
shown in Fig. |^. The decay through the virtual photon or Z (Fig. |^(a)) contributes only 
to the channel W~^W~, the decay through the virtual Higgs (Fig. |^(b)) is shown only for 
completeness. It contributes neither to y nor to 0 decays because of its quantum numbers 
(this only applies to a standard model Higgs, of course). 

First order QCD corrections receive contributions from the diagrams shown in Fig. |^. 
Their sum is infrared hnite. Real gluon emission is forbidden by color conservation. Compar¬ 
ing this result with the calculation of QED corrections, which is practically the same in this 
order, we could also argue that the coupling of real gluons (photons) to a color (electrically) 
neutral state vanishes in the static limit, thus providing an infrared hnite answer. 

However, as expected, all corrections exhibit the Coulomb singularity, e.g. 

r(y ^ 77) = F(®°™)(y ^ yy) + Y 



which is universal, proportional to (Cir = 4/3 is a color factor), and which originates from 
box- and from s-channel vertex correction diagrams (Figs. |^(a) and (d)). This divergence 
actually represents part of the Bethe-Salpeter wavefunction and must be dropped since it is 
already included in the factor ^^(O). Furthermore, since we are calculating ratios of decay 
widths, the singularities (formally) cancel anyway. The K-factors for the ratios R, and their 
nontrivial parts 6k, which are dehned through 


^ab — 


corrected ' (^“0 


Born 


= 1 + 


271 


6k^ 


( 5 ) 


will thus be free from Coulomb singularities. The corrected rate can then be obtained from 




Born 


■ ^ 77) ■ (l + 


asCp /TT^ 
271 V 2 


10 + 


( 6 ) 
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FIG. 1. Tree level generic Feynman diagrams 


s-channel corrections 




Z“(Y) 


(b) 


t-channel corrections 



(c) 


(d) 




(e) 


(f) 


FIG. 2. Diagrams contributing to first order QCD corrections 
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At this point a comment on the regularization the Coulomb singularity is appropriate. 
Two different procedures are possible: one may either start with nonvanishing (3 and consider 
the limit /3 —>■ 0 in the end, which obviously requires signiffcantly more effort during the 
calculation than really needed, in particular for the box diagram. Alternatively one may set 
/? = 0 from the outset and employ a nonvanishing gluon mass A. This second procedure has 
several advantages: the Coulomb singularity and the infrared divergences are regularized 
in one step and the special kinematical situation facilitates the calculation (especially of 
the box diagram) significantly. To connect the two approaches the vertex correction can be 
investigated. This leads to the substitution rule 

m ^ TT 
A 4|ff|’ 

Before the results of the calculations can be presented, the notation must be fixed. Vector 
and axial vector coupling constants are abreviated with the help of 

V = 2{Isl + hn) — 4(5 sin^ , a = 2{I3l — , ?/ = 2 sin 2Qw 

where Q denotes the quark charge divided by the proton charge, I 3 (the 3-component of) the 
weak isospin and the weak mixing angle. The results are applicable to fourth generation 
quarkonia as well as to more unconventional quarks (for example an isosinglet in Eq models) 
as long as the couplings to the Higgs boson coincide with those of the standard model. 

All masses are measured in units of the quarkonium mass M = 2mq, whence px = 


III. RESULTS 


A. Decays into ZZ or 77 


Decays into two Z bosons are possible if the mass M of the quarkonium is larger than 
twice the Z mass and proceed via diagrams |^(a,c,d). -0 77 is forbidden by C conservation, 

whereas 5 ’ ZZ is allowed through the parity violating coupling. T{ 7 ] ^ 77 ) can be 
obtained from r (?7 ^ ZZ) by replacing a —> 0, v ^ Qy and taking the limit pz 0, which 
leads to the result given above (|[). The lowest order predictions for the ratios 


{Qy)^ (1 - 2pz)^ 

pP ^ 2 (an)^ ( 3 ^ 

^{Qy)^pz{l-2pz)^ 


( 7 ) 

( 8 ) 


where (3 = — 4M|/M2 in this case, are well known. To include the QCD corrections, 

they have to be multplied by the K-factors, which contain as their nontrivial parts: 




(3^ 


TT^ 

y + 


1 + 3/32 - sy -f / 3 « 

y 


ln(l + y) -8 


2 + 5/32 - (3^ 

y(i+y) 


B{pz) 
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l + /9^ 


C{(3) + 


8a2(l + /32) 


(3^ {v^ + a?)(3‘^ 




Skt, = 2 - ^ _4 

+ P^)C{P) 


.{1 + P^)C{P) 


6 + 7/32 + 11/34 + 2/3® 

P^OTW) 


B{pz) = \Jpz{l - pz) arcsin 


2/3 4 


+ Lio - 


Li2(a;) denotes the dilogarithm; Li2(a;) = — fn . 


+23? Li 


1-/3 
1 + /3 
1 + /3 





The limiting case of a large quarkonium mass as well as shorter approximate formulae 
are presented in section IV. 

The two K-factors are shown in Fig. |^, using 

127r 

The leading apparent singularities of K'^z and K^zz at threshold (M = 2Mz or /? = 0) 
proportional to l//3^ and l//3^ respectively cancel. Hence K'^z diverges oc l//3^, K^zz remains 
hnite: 


Sk^z 


5k%z 


29 TT 
v/Stt 16 


TT^ 1 

dSvr 


x/Stt 8a^ 

3 v'^ + a? 

2 


/3 ^ 0 


5^3 


— (3 

2 


0 


( 11 ) 

( 12 ) 


The D-wave phase space oc for 'ip ^ ZZ which is present in Born approximation is 
thus modihed to a behavior oc j3'^ characteristic for P-wave phase space. 


B. Decays into 

The decay to a Z boson plus a photon is possible if the quarkonium mass exceeds the Z 
mass. The t- and u-channel quark exchange diagrams (c) and (d) contribute. The normalized 
decay rates are given by 


Rl, = 2 


{Qy) 


r/? 


^ 2 m-(^) 

3 {Qyf (1 - 

(with j3 = 1 — and the K-factors by 

.2 1 n A a or2 ^ 


(13) 

(14) 


Ap _ g 1 Ttt 1 - /3 4/3 - 3/3 

l - 2;5 6 13 ( l - 2/?)2 




.3-2/? . 2 

2-::— arcsm ^pz -|- 2 


ln(2/?) 

1-/9 


5kt,= 


2(3 


+ 


/9 

TT^ 16 - 24/? - 5/32 ^ 7^3 


/9 


Li2(l-2/?) 


(15) 


-4 


1-2(3 ' 6 
8 -5(3-(3^ 


/?2(2-/?) 


^8-39/? + 63/?2-32/?3 + 2/?^ 
^- P{ 2 -P)(l- 2 pr - 


B{P2 


/?2(2-/?) 

8 - 15/? + 2/?22/33 

/92(2-/3) 


^ 8 - 13/3 + /32 + 2/33 , 2 ^ 

Li2(l-2/3) 


(16) 


with B{pz) = \JpziX — Pz) arcsin ^/p^ as before. 
These two functions are shown in Fig. 0. 
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FIG. 4. Kz^ and 5kz^ 

Close to threshold /3 approaches zero and diverges, in contrast to Sk^. Specihcally: 

(17) 

Ski, ^ “ 7 ■ ^ 

This behavior resnlts from the combination of the singnlarities of the glnon and the qnark 
propagator. In fact, the virtnal quark is close to its mass shell in this limit. Very close to 
threshold intermediate bound states should be taken into account since an alternative way 
to describe the decay is through the chain 1 — Z) with mixing between I’*"’'' 

and the Z. A similar problem also appears in yiif [Q. For (3=1/2 the correction is 
regular. 
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C. The decays into ZH and 


As stated in the introduction the decay rj ^ ZH could signal the presence of the Higgs 
boson and the new quark simultaneously. In Born approximation one obtains for the ratios 


R: 


R 


V _ 
ZH — 


a? j3‘^ 

{QvYpI 


2 v^ (3 

KQvY Pz{l - Pz - PhY 


X 



— pz — Ph 


1 — ^Pz 
^ — Pz 


2 



2 

- Pz + Ph- 


- pH 

- Pz 


2 


( 19 ) 


( 20 ) 


with (3 = \J3 -P p\ + — 2pz — 2pH — 2pzPH- In principle, all three diagrams |^(a,c,d) 

could contribute to both decays. However, the contribution from (c) and (d) vanishes for 
p ^ ZH (this is no longer true for the QCD corrections), whereas R^h is the sum of three 
terms: 


jdP 

^ZH,t 


T3p _ 
^ZH,s — 


jdP 

^ZH,st 


2(3 (3^ + Qpz 


{QvY ^pz {i- pz- phY 

PP^ + 12pz 


t-channel, (c) + (d) 


{QvY 3 (1 - PzY 

8(31 + Pz - Ph . 


{QyY 3 1- pz + pH 


s-channel, (a) 

interference term 


( 21 ) 


For large quarkonium masses, RPzh rises proportional to l/p|, in contrast to all other 
modes listed up to now, which increase proportional to 1 /pz or approach a constant value. 
The different behavior can be understood from the Goldstone boson equivalence theorem 
0 which identihes the amplitudes for scattering of longitudinal gauge bosons {Z or W) at 
high energies with those for the appropriate Goldstone bosons. The latter must be taken as 
pseudoscalar particles. The decay p —> GzH {Gz denotes the Goldstone boson belonging to 
Z) is allowed, i.e. the Z emitted can be longitudinal, with its coupling proportional to M/Mz- 
This explains one factor Ijpz = The second one originates from the coupling of 

the Higgs to fermions. In contrast the decay Y GzH is forbidden by GP conservation, 
and only the transverse part of Z with a coupling proportional to Pweak remains. This leads 
to a rate oc 1 /pz- 

The behavior of the other rates can be obtained through similar arguments: p 
GzGz^GzZ or Gzl are forbidden, hence RJ^z RJ^^y approach constant values. The 
decays Y GzZ (but not GzGz) and Gz'y are allowed, these ratios increase cx 1/pz- 

The full expressions for the K-factors for Rzh can be found in the appendix. The 
numerical results are shown in Fig. |^. The corrections are large, at least about 15% for R^zh 
and 30% for Rzhj but of course they do not drastically change the conclusions that can be 
drawn from the Born results, at least if we are not too close to threshold. 













1.5 1.75 2 2.25 2.5 2.75 3 3.25 3.5 4 



FIG. 5. Kzh and 5kzH for Higgs masses Mh = 70,100,150,200GeH (solid, dashed, dotted, 
dash-dotted). 

In the limit 0, keeping pn = fixed and also keeping a factor Ijpz = 

{ygn/^Y in the numerator, the well-known result for the decay V’ 'jH is recovered 0]: 


ph _ 




(l - Ph) 


( 22 ) 


„ l-3)3 Tr'‘ „ 2-/3" ' , 4-13/3+ 7)3^+ 8)3" , 

^ 2 ^ /3(1-2/3)2 


1 /7r2 

+ -(-Lia 1 - 2 

/3V6 


+ 4 ^ ^ + Li 2 (l - 2/3) - 2arcsin^ 




(23) 


where /3 = (1 — Mjj/M"^) and i3(p) = \lp{3 — p) arcsin In the same limit one obtains 

1-4/3 7r2 8(2-/3) , ^^2 , 1-/3 






'1-2/3 


(1 - /3)/32 
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+2^-^(y + Li2(l - 2^) - 2arcsin2 (24) 

The limiting behavior of is well known: 

5A:^^ = -4-121n2-^ {(3^1) (25) 

StT 1 TT^ 

Sktn = - + 8 - y (0^0) (26) 

For Sk^Hj 0^6 obtains in the corresponding limit 

Sk^H = -2 - 8 In 2 {(3^1) (27) 

Skin = -y ^ + ® - y ^ P8) 


D. The decays into W~^W 


Although the decays into W^W~ seem very similar to the decays into two Z, important 
differences arise: two distinctively different isospin assignments must be considered: case 1 
for an isosinglet quark and case 2 for a standard model quark with I/ 3 I = 1/2 (isodoublet). 
In case 2 the propagation of the isospin partner of the quark forming the bound state enters 
the diagram. Its mass is generically denoted by my. For an isosinglet quark (case 1), only 
diagram (a) is possible, and we obtain 


f?k! _ n 

■‘^ww ~ 

ph _ ^ 

^ww — 


f3^ 


Q 2 48(1 - pzY ^ ^ 

In case 2, diagrams (a) and (c) contribute, and we find 
1 

W ^__ 

(Qsin 0 )^ 8(1 + 4pT - 4pty)2 


Mz 


M, 


w 


pb _ 

^ww — 


(3^ 


(Qsin 0 )^ 768p^ 


( 2 / 3 ) 


2(l + 20p^y + 12p^)(l-4cpz)' 


il-Pzf 


,rT {l-Acpz){^ + Qpw) , 2 


2 — pw 


(1 + 4pT — 4piy)2 


(29) 


(30) 


(31) 


where 

^ . 9 ^ 

= Y 1 — 4pvv , Pt = , c = 2/3 sin 0^ 

The dependence on the sign of which is in particular relevant for the interference term 
is displayed explicitly^], so our result applies to up- as well as to down-type quarks. In the 


iThe sign of this term is in conflict with the formula given in for d-type quarks. 
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limit of a large quark mass, whereas approaches a constant. This 

again is a consequence of the equivalence theorem, which tells us that in the "0 decay both W 
can be longitudinal. In contrast both rj GwGw and r] —> GwW are forbidden. Obviously 
the WW mode would dominate 0 decays. 

For case 1 (isosinglet) the K-factor is simply constant: 

^^t/w = 2 — (32) 

The K-factors for case 2 (standard model) are again quite lengthy and are listed in the 
appendix. Formally = Sk^z if we set ttit = m, v = a. Fig. shows the two functions 

for three different values of m^. 



FIG. 6. 5Kww and Kww for down-type quarks and rriT — m = 0, 50,80GeV (solid, dashed, 
dotted). 
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IV. LARGE QUARKONIUM MASS AND APPROXIMATE FORMULA: 


The K-factors approach quickly their asymptotic behavior for (3 towards zero or 1. Since 
additional quarks are presumably rather heavy, it often will be sufficient to take the limiting 
values px —^ 0. They are given in table |I[ In the case of WW a possible mass splitting 
between the constituent quark and its isospin partner is neglected. 

Up to obvious coupling constants, the ratio between ^-nd approaches 1/2 in 
this limit, a consequence of the statistical factor for identical particles in the ZZ case. The 
K-factors therefore approach the same value. Also and become equal, and the 
asymptotic value of remains unaffected by the QCD corrections. The difference between 
only arises from their different couplings to left- and right-handed quarks 
and the appearance of rriT in the WW case. It should be mentioned that - for a 

down-type quark - approaches its asymptotic value very slowly, so in this case the full result 
given in the appendix should be used. 

The following approximations are valid for the full kinematical range: 

5klz ~ -5.84 - 4.36pz - 47.73p| (33) 

/o 

Sktz ~ - 13.66 + 3.09pz - 35.41p| (34) 

1 - Apz 

5kl^ ^ -0.642pz - 0.473p| (35) 

5kt,^-^-^^^= + Rl (36) 

3 Vl - pz 

They reproduce the exact result within an error of less than 1% even close to threshold. 


TABLE I. Asymptotic values of the (5A:-factors for px —> 0. 
Note that in order to obtain the full corrections, these numbers 
have to be multiplied by as{M) ■ CfI2tx. 


K-factor 

O 

T 

numerical value 

Sk^z 

-^(vr2-81n2) 

-5.84^^ 


0 

0.0 

6k^H 

-2 - 81n2 

-7.55 


—TT^ -|- 81n2 

-4.32 

Sk'zz 

_^ + 2-41n2 

-8.17 

“h 

_^ + 2-41n2 

-8.17 

5k^H 

-4- 121n2 

-14.79 

^kf^H 

-4- 121n2 

-14.79 

^kww 

2(1 + 2/3 )-^(1 + 4 / 3 ) 

4.93'’’'’ 



-10.80‘’’’’ 


^down-type standard model quark assumed 
'’equal masses for the isospin partners assumed 
‘’up-type standard model quark assumed 
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V. SUMMARY 


The complete evaluation of the QCD corrections to annihilation decays of superheavy 
S-state quarkonia into elektroweak (gauge) bosons has been presented. Nearly all the correc¬ 
tions are negative, with -0 —>■ WW (d-type quark) being the only exception. Some of them 
are sizeable even in the high mass regime, where all mass scales can be neglected and the 
K-factors become simple constants. In particular the corrections to "0 ^ ZII and yif are 
of the order of 30% in this limit, whereas the others are all below 18%. The dominance of 
the mode rj ZH among the 77 -decays for a certain range of Higgs and quarkonium masses 
is not changed by strong radiative corrections. 


APPENDIX: 


1. Generalties 


Before discussing the full analytic results for Kzh and Kww, the following functions 
have to be introduced: 


Qoi — ( 


2 2 

m — 




fim — mj'\2 


M 


H 


X arctan 


1 _ f m-mr \ " 

\ Mh J 


A / m-\-mT A " 
Mh ) 

Qo 2 = Qoi {Mh —>■ Mz 


2 2 2 2 
^ Trim — q . — q 

Qi 2 = -^- hi-^— 




rrin 


In^ 


Qi = rn^Colrn^, M|, 0, m, rriT) 


Q2 = rri^Co{rri^, M^, 0 , m, itit) 


Qo = 4m^, Mjj, g^, M|; A, m, m, tut) + 


m 


2 A2 TT^ 


2 {rri^ — g2) ^ m? 


(^In 


+ 


A is the gluon mass to regulate the infrared and Coulomb divergences as mentioned in section 
II, V denotes the quark velocity and the scalar three and four point one loop integrals Cq 
and Dq are defined by m 


Do{pI, (pi -P2)^pL {P2-P3f,pl, (P3 -Pi)^;nro,mi,m2,m3) 

_ _ 1 _ 

J [/2 - ml] [{I + piy - mj] [(/ -h ^ 2 )^ - ml] [{I + p^y - mj] 


/ (^ 
- 

i 


dH 


[/2 _ ml] [(/ + piy - mf] [(/ + P 2 y - ml] 

In the case of Kzh, mx = m, which makes the expressions more compact, and for Kww 
both Mh and Mz must be replaced by Mw- 
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When taking the limit A —0 the four point function Dq develops an infrared and a 
Coulomb singularity which are exactly of the displayed separately in the dehnition of Qq, so 
the latter is hnite. When calculating Qo, it is not necessary to work out the full expression for 
Do, because use can be made of the special kinematical situation which results in relations 
between the propagators appearing in the denominator: 

2A2 = [(/ + - m^] + [(/ - - m^] - 2[/2 - A^] 

2/2 = [(/ + - rri^] + [(/ - - m^] 

where Pq is the momentum carried by the incoming quark and thus, in our approximation, 
—Pq the momentum of the antiquark. These identities can be used to rewrite Dq: 


Da = 


Mjj, {pq - pzf] A, m, mx) - Co{m^, {pq - pzf] 0, m, mx) 
+Co(m^ M|, {pq - pzf; A, m, mx) - Co{m^, M|, {pq - pz)‘^; 0, m, mx)) 


2A2 


and the calculation is reduced to the problem of expanding three point functions, with the 
result: 


Qo{ZH) — — d{pz) + d{pH] 


d{Px) - 


vApx + z-2f3 


with 

z 


(Al) 


z — 2(3 2 — z — 2px 

Apx - 1 + 2(3 jz 


2(3 8„, , 2{l-z) 

In — - -B{px) + - -^ 

z z z 


In ■ 


4/9(1 - Ijz) 


4px + ^ ~ 2 — 2/5 + 4/5/z z + 2/5 + 4px ~ 2 

1 r 2(1 - 2 yiz) 


In 


Apx + ^ — 2 + 2(3 
4/9(1 - l/z) 


In 


z(yZ + 2(3 + Apx — 2)1 

A(3lz 


z-2(3^Apx-2yz-2(3-VApx-2^ A(3lz 2(3 + 2-Apx-z 


z + 2 — Apx — 2(3 


In 


4/9 


^ — 2(1 — pz — Ph) 


Qo{WW) = --Jn^ 


arctan 


z-2 + Apx + 2p 2p + 2-Apx - zi_ 

{2 - z){2e - z) + (3{At‘^ - zy 


y/W 


(2/2 


arctan 


/9(4/2 


^4/2 — 2/2 (2/2 — z){(i + (3) 


V4/2 — z‘^{2D — z — (3{2 — z)) 


z = e + (3\ t<i + 


Ml 


w 


m 


(A2) 


To demonstrate the hniteness of Qq, a third relation can be used, which allows to ex- 
plicitely split off the divergent part of Do(A —0): 

f + 2l{pq - Pz) - [(/ + Pq- pz)"^ - rny\ =mT- 


resulting in 

Do(A ^ 0) = 


m\ — g2 


Co(m^, 4m^, A —0, m, m) 


■dH 


/2 + 2l{pq-pz) 


m"^ — J P [(/ + Pq)^ — m2] [(/ — pq)"^ — m2] [(/ + pq — px)"^ — rnx 


2m2(m|. 


■ f In —- + 


TT 


2\v\ 


— 2) + hnite integral 
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This equation, however, is not of much practical use because the calculation of the remaining 
hnite integral by standard methods would re-introduce infrared divergencies when expanding 
it into two- and three-point functions. 


2. K-factors for the ZH mode 

With these ingredients, the missing K-factors can be calculated. First the results for 
Kzh'- 




3 - ipz - 4 :Ph . f Pz - pH 

l-2pz-2pH 2^1 (3^ 


+ 


1 — 2 ,pz — 2 pH 


2Qi2 + 4— 


+4|q„. - ,A3) 

The only form allowed for the amplitude for p —> ZH is proportional to ez-Pqq- Hence the 
radiative corrections can trivially be written as a multiple of the Born result and the K-factor 
assumes a fairly simple form. However, the same consideration does not apply to 'll! ^ ZH 
where two quite different and more complicated structures are present for the amplitudes on 
the tree level. Therefore the corrections cannot be written as a single ‘compact’ K-factor. 
However, the explicit calculation shows that every amplitude corresponding to a given graph 
of Fig. 1^ can be split into two parts that are each proportional to one of the Born amplitudes. 
Hence, after multiplication with the latter and after spin summation, the three expressions 
m are recovered and ‘partial K-factors’ can be read off (with z = 2(1 — pz — Ph))'- 


6 ks = 


5k St = 


1 + Pz ~ 2zp^ 
Pz(l - z) 


^ _ (1 - pz){2 - 3z) 
2 pzz(l - z) 

pz ^ ’ 


2(1 


2 2 z{\-z 

r3pz -I- pH 


Qi2 + 2Fi{pz, Ph) 


Qi2 + Fi{pz, Ph) + F2{pz, Pir)d 


(A4) 

PzFsIpz, Ph) 


12 


6pz(l + Pz — Pir)2 
4(3pz — 1 + Ph) + ^(1 + Pz — Ph) 


Qi 2 + “ l)Qoi — PzQo 2 + 8 pz{z -|- 3pz)Qo 


^ ^ Spz + z{l + Pz - Ph) - 2(3“^ ^ 

VI H-]-V2 


5kt = 


A 42 2 ^ C)^ 

-7]-^Qi2 + 2F2(pz, Ph) + ^2 4c— ^sipz, Ph) 

z{l - z) + Qpz 

n2 a ~ Qoi + ‘^PzQ2 + (t ~ l)Qi + 2(2/3^ -|- 12p^ -|- z)Qq 
f3^ + Qpz ^ 4 


z vr 

~z~^ 

z 


(AS) 


(A6) 


where 

Fi{pz,Ph) = 


Pz 


((1 + Ph — Pz)Qoi — 2pzQo2 + 


Qpz — ^ — Ph 


PzP 2 V- ■ -- ■ 2 

+ [{Ph — PzY — l]<5i + 2pz(l + Ph — Pz)Q‘^ 


Q 


12 


Ph) = 1- %iQoi + + 4 






I- pH + 


Ph(1 — pH + Pz) 




Qi 
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(3“^ + Apl; - Apz + l{4:pz+ 2 - z)+ 2pH [i , ^ , ^‘^PzPh 

H-- Qi 2 - 1 + 4p^ + 


^2 L ' ' (32 

Fsipz, Ph) = ^ [12(pz - 1 - Ph) - 4^]Qo2 + 72 ^ _^^ Qi2 

/3^(64p2 - 20zpz + z^) + zpz{96pl “ 216pz - 12 + 242:p2) ^ 

H- r^TTZ -^-V12 


4p|(l-z) 




+2[l-(6 - z) - 12z + 24(1 - pz)]g 2 
Pz 


Q 2 


The corrected ratio is then given by 


F'zh — {J^h) 3 {^ksRzH,s + ^^stR'zH,st + ^^tRzH,t) 

It shonld be stressed, however, that the notation does not imply that for example 6 ks is 
the correction induced by the s-channel diagrams £g. 0(a,b) only. In fact, it contains part 
of the t-channel contribution and must be considered in the limit Mz —> 0 to obtain the 
correct result for although on the Born level only the t-channel diagrams are possible 
in this limit. 


3. K-factors for the WW mode 


The last missing K-factors are those for the decay modes "0 —> WW and p —> WW. 
They read 


„r, (3^1 - (3^ - 2 R) - A(3H^ - (3H^ - R + /3^ ^2w2^ 

^klw = —3-2 -V-^-In' ' 




P^ 




+ \n{f 


+2 I 32 oA _ 

_(2 -R- + 2'-4P^Q2 + 2 ^^Q, 


P^ 

For there are again three partial K-factors: 


P^ 


(AT) 


TT 


5fco = 2- 


5h = 


{p^ + R){R + R - 2I3R TT^ R{z^i3^j^aR^ + 7R + RH^ + R) ^ 


P\7 + PR 


PR7 + PR{P^ + R) 


HR 


+ 

+2 


H^{P^ - 4fg -t ^ QtP - 11/g + 2 p^ + 2 tp^ + 2 RH - 


PRPHR) PH7 + P^) 

Af3^ - 4/32 -(3H3t- 2tp^ - tH + 4f2 - R(3^ - R(3^ - R 


Q 


12 


go2 


(A8) 


PR7 + PR 

Ap^+R{U-6t - ISR + 3R)+P'^{U+ 6t - 2AR - 2R + 2R)-p^{3 -2t- SR) 


PR7 + PR 

14/32 - 9/3^ + 4/36 - 32(14 + 8/32 - 7/3^ - 5R - 2pH^ + R) ^ ^ 


Q 2 


PR7 + PR 


iP^ + R)Qo 
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6 k.f = - 


TT^ (/32 + t2)(13 - 13/32 + 332 + 2/3^ - 32/32) 

T 


+ 


2/32(1-/32)(13-3/32) 

32(130/32 - 142/3^ + 29/36 + 2632 - 832/32 + 832/34 + 634 + 334/32) 

2/32(1 - /32)(13 - 3/32)(/32 + 32) 
r/32(13 + 143 - 232 + 534) - /34(27 + 163 - II32) + 2/36(8 - 3) - 1332 - 334 


■ ln(t) 


232 


Q 


2/32(1-/ 32 )( 13 -3/32) 

73 + 34 - 10/32 - 83/32 ^ 4^2^2 ^ 9^4 ^ ^^4 


^2+/32J^12 , /32(13 - 3/32) 

r/32(26 + 143 - 8232 + 1633 - 1734 - 336) - ( 3^{58 + 163 - 37 f + 23^ + 34) 


Qo 2 


+ 


2/32(1 -/32)(13 -3/32) 

/36(27 + 23 + 32) - /33 + 2632 - 1433 + I634 - 336 - t 


+ 


Q 2 


2/32(1-/32)(13-3/32) 2/32 

/32(26 + 1032 + 534 + 36) - /34(51 - 1732) + /36(25 - 3 t ^) - 2/3® - 2632 - 334 ^ ^6 


2/32(1 -/32)(13-3/32) 

X (/3^ + t^)Qo 

The corrected rate can be obtained analogously to R^h- 
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